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FOREWORD 


This  report  describes  work  aimed  at  investigating  the  behavior  of  composites  under 
multiaxial  loading.  Work  on  the  project  began  at  the  University  of  Houston.  Experiments 
to  determine  local  variations  of  deformation  and  local  mechanicalproperties  have  beed 
developed  and  demonstrated.  In  August  2000,  the  Pi’s  accepted  a  faculty  position  and 
moved  to  the  University  of  Texas  at  Austin.  This  resulted  in  an  interruption  of  the 
project;  the  project  was  terminated  at  the  University  of  Houston,  and  then  reissued  to  the 
PI  at  the  University  of  Texas  at  Austin.  This  process  took  about  eight  months  to 
accomplish.  In  this  report,  the  work  done  at  the  University  of  Houston  is  described. 
Further  work  on  applying  these  methods  to  the  examination  of  the  deformation  and 
failure  of  composites  is  currently  underway  at  the  University  of  Texas. 

STATEMENT  OF  THE  PROBLEM  STUDIED 

We  believe  that  through  an  understanding  of  the  micromechanical  deformation  and  load 
transfer,  design  with  composites  can  be  approached  by  coupling  the  microstructural 
design  and  macrostructural  design.  This  approach  will  eliminate  the  need  for  ad  hoc 
strengths  and  provide  for  design  using  the  fundamental  material  properties  of  the 
constituents.  In  order  to  develop  this  methodology,  it  is  necessary  to  understand  the 
deformation  and  failure  behavior  of  composites  at  the  level  of  the  heterogeneity.  This  is 
the  objective  of  the  present  work.  We  are  examining  a  multiaxial  loading  condition  as 
described  in  Figure  1 .  While  most  experiments  focus  on  measurements  at  the  global 
level,  we  combine  these  global  measurements  with  measurements  of  the  local  variations 
of  the  deformation  through  a  digital  image  correlation  technique.  Local  mechanical 
properties  are  also  necessary  to  implement  these  ideas.  With  a  view  towards  establishing 
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this  experimental  strategy,  we  have  two  experimental  schemes  that  are  currently  under 
development.  First  is  a  digital  image  correlation  method  for  the  measurement  of 
displacements  with  a  resolution  of  20  microstrains  over  spatial  domains  that  are  on  the 
order  of  a  few  hundred  microns.  The  second  is  a  technique  for  evaluating  the  properties 
of  the  constituent  materials  under  multiaxial  stress  states  similar  to  that  experienced  by 


Figure  1.  Macrostnictural  geometry  and  loading  configuration.  Multiaxial  loading  and  out-of-plane 
defomiations  can  be  varied  by  varying  the  parameters  in  tliis  configuration. 
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SUMMARY  OF  IMPORTANT  RESULTS 

We  describe  briefly  the  results  obtained  using  the  digital  image  correlation  method  and 
the  confined  compression  scheme  for  multiaxial  characterization  of  materials. 


DIGITAL  IMAGE  CORRELATION 

This  method  has  been  used  by  a  number  of  investigators  (for  example,  Peters  and 
Ranson,  1980,  Sutton  et  al,  1991,  Vendroux  and  Knauss,  1998)  The  principle  of  the 
method  is  illustrated  in  Figure  2,  where  the  basic  analytical  element,  the  subset  size  is 
defined  as  a  region  of  ±A  pixels  around  the  reference  point  xo.  A  comparison  of  small 
subsets,  centered  at  those  reference  points,  yields  local  strain  and  displacement 
measurements  at  those  pre-defmed  points.  This  is  discussed  in  the  following  for  one¬ 
dimensional  deformations,  the  formulas  are  generally  applicable  for  the  three- 
dimensional  case.  Let  y~f(x)  represent  the  deformation  and  let  A(x)  and  Biy)  define  the 
intensity  variations  of  the  images,  before  and  after  deformation,  respectively.  We  first 
define  an  array  of  points  on  the  undeformed  image  where  the  displacements  and  strains 
will  be  computed.  Then  we  define  a  subset  S  of  the  undeformed  image  around  each  of 
those  reference  points  xq.  To  that  subset  is  applied  a  function  y(x)  that  is  supposed  to 
‘model’  the  deformation,  such  that 


/(x)  =  x  +  m(x„)  + 


(X-Xa) 


(1) 


The  numerical  intensity  A(x)  of  the  undeformed  image  at  points  x  of  5  are  compared 
to  the  intensity  B{f{x))  via  a  correlation  coefficient  C.  The  set  of  displacements  m(xo)  and 
strains  [clu/dx]xo  that  minimize  the  correlation  coefficient  becomes  the  solution  for  the 
displacements  and  strains.  The  cross  following  correlation  coefficients  has  been  widely 
used; 


iAix)B{f(x))dx 

S  _ 


j  A{x)'dx\  B(f  (x)y  dx 


1/2 


with  S  such  that  x  e  5"  o  (Xg  -  A)  <  x  <  (x^  +  A) 


(2) 
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Profile  A(x) 


One  p4Cture 
element  (Pixel) 


Figure  2:  Principle  of  tlie  Digital  Image  Correlation  method  -  Detail  of  an  image  showing  the 
decomposition  in  pixels  -  [|:  Picture  of  the  undefonned  microstructure  -  Spatial  gray  scale  intensity 
variations  A(x)  over  a  horizontal  line  of  pixels  of  the  undeformed  image  1  -  Picture  of  the  deformed 
microstructure  under  a  tensile  load  F  (the  deformation  is  exaggerated  for  effect)  -  The  dashed  line  shows 
the  intensity  variations  B{x)  over  the  same  line  as  in  2.  The  solid  line  represents  the  corresponding 
interpolated  intensity  profile  -  An  assumed  strain  is  applied  to  the  subset  S  of  1  through  the  function 
/l(x).  Its  image  through  B{f{x))  is  compared  to  the  undeformed  intensity  profile  A{x)  for  various  values  of 
the  strain y;(x),y^(x)  etc. 

In  order  to  get  subpixel  accuracy  the  subsets  from  the  deformed  images  must  be 
interpolated.  The  interpolation  method  has  been  shown  to  influence  the  accuracy  of  the 
results.  Although  increasing  the  subset  size  generally  improves  the  quality  of  the 
comparison,  the  effect  of  averaging  the  strains  over  the  subset  and  an  increased  CPU  time 


must  also  be  considered.  The  successive  comparison  of  the  subsets  involves  an  iterative 
method.  The  first  method  used  was  the  trial  and  error  method,  where  the  subsets  are 
compared  over  each  point  of  a  given  area,  sequentially  reducing  the  size  and  the  step  of 
the  search  around  the  point  of  minimum  correlation  coefficient  until  the  target  accuracy 
was  reached.  This  method  is  time  consuming  and  later  improvements  yielded  to  the  use 
of  a  gradient  method  (Newton-Raphson  scheme)  which  considers  the  gradient  of  the 
correlation  coefficient  in  order  to  accelerate  convergence.  Unfortunately  this  method  can 
stall  at  local  minima  and  hence  miss  the  solution.  The  correlation  coefficient  must  also  be 
considered.  Cross  correlation  coefficient  has  been  commonly  used  but  other  coefficients, 
such  as  a  least  square  coefficient  gave  better  results.  This  coefficient  can  be  written  as 
follows 

(-=5 -  (3) 

\A{x)'dx 

s 

The  accuracy  of  determining  the  displacements  by  a  the  digital  image  correlation 
method  was  examined  analytically  and  experimentally;  these  are  fully  documented  in  a 
manuscript  under  preparation  (Fossier,  Ravi-Chandar  and  White,  2001).  Here,  we  show 
the  capability  with  two  figures,  one  illustrating  the  method  in  a  ceramic  and  the  other  in  a 
polymer.  Extensions  to  composites  are  quite  straightforward.  Figure  3  shows  the 
application  of  the  method  to  an  magnesium  aluminate  spinel  sample;  a  microtensile 
sample  was  loaded  and  viewed  under  a  microscope;  a  surface  coating  was  provided  to 
enhance  the  DIG  method.  The  elastic  behavior  under  small  loads  (within  a  single  grain) 
and  the  Poisson’s  ratio  (also  within  a  single  grain)  were  measured.  Clearly  the  DIG 
method  provides  measurements  of  mechanical  properties  at  the  scale  of  a  grain.  Figure  4 
shows  strain  contours  obtained  for  a  polycarbonate  sample  in  a  three-point  bend  loading 
over  a  field  of  view  of  800  pm  x  800  pm;  while  the  global  measurement  of  strain  is 
uniform,  the  local  measurement  indicates  significant  nonuniformity.  The  nonuniformity 
of  the  strain  distribution  is  evident  in  this  figure  and  is  a  result  of  structural 
nonuniformity.  These  examples  were  a  part  of  our  initial  studies  to  make  sure  that  the 
image  correlation  method  was  appropriate  for  the  task  of  determining  local  strain 
distributions  in  a  composite  specimen. 
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Figxire  4.  Strain  distribution  in  a  polycarbonate  specimen  under  three  point  bending.  Contours  indicate 
constant  strain  levels.  Obtained  using  digital  image  correlation 

MULTIAXIAL  MECHANICAL  PROPERITES 

In  the  composite,  the  constituents  are  under  multiaxial  loading  conditions  even  when  the 
global  loading  is  a  simple  uniaxial  loading.  In  order  to  determine  the  actual  local 
mechanical  behavior,  we  have  developed  the  confining  cylinder  experiment  shown  in 
Figure  5.  The  specimen  is  a  cylindrical  rod  of  diameter  2a  and  length  L.  It  is  inserted  into 
the  cavity  of  a  hollow  cylinder  of  internal  diameter  2a  and  outer  diameter  2b\  ideally  the 
clearance  between  the  specimen  and  the  confining  cylinder  should  be  zero.  It  is  assumed 
that  the  deformation  in  the  cylindrical  specimen  is  homogeneous.  The  material  of  the 
cylinder  is  chosen  such  that  it  provides  the  necessary  elastic  restraint  to  the  specimen  and 
deforms  elastically  during  the  entire  deformation  of  the  specimen;  the  cylinder  must  also 
retain  sufficient  compliance  so  as  to  enable  measurement  of  its  deformation  during  the 
test.  The  cylindrical  specimen  is  strained  by  inserting  a  plunger  into  the  cavity  and 
displacing  it  by  a  loading  frame  to  generate  a  uniform  compressive  strain  : 

e^(r,0,z)^ea.  (4) 
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The  cylindrical  polar  coordinate  system  used  is  shown  in  Figure  5.  The  corresponding 
axial  component  of  stress,  a  ^{r,6,z)  =  cr„,  is  uniform  in  the  entire  specimen  and  can  be 

measured  easily  using  the  load  cell  in  any  standard  tensile  testing  machine. 


Figure  5.  Configuration  of  the  confined  compression  apparatus  for  multiaxial  characterization. 

The  confining  cylinder  provides  a  restraint  against  the  radial  expansion  of  the  specimen, 
thereby  generating  a  radial  pressure  on  the  inner  surface  of  the  confining  cylinder.  The 
continuity  condition  on  the  lateral  surfaces  of  the  confining  cylinder  and  specimen  may 

be  written  as  follows:  (a)  =  (a)  if  and  only  if  (a)  =  (a)  >  0  and  (a)  =  0 

otherwise*.  The  main  idea  of  this  type  of  constraint  is  to  provide  stability  to  the 
deformation;  in  the  absence  of  this  elastic  restraint,  barreling  deformations  can  and  do 
develop.  Note  that  the  radial  pressure  may  also  be  provided  by  performing  the  test  in  a 
liquid  under  high  pressure;  there  are  a  number  of  experiments  in  the  literature  that  deal 
with  the  response  of  materials  under  hydrostatic  pressure,  beginning  with  the  work  of 
Bridgman  (1964).  However,  hydrostatic  pressure  from  a  liquid  medium  does  not  stabilize 
the  deformation;  in  other  words,  the  hydrostatic  pressure  acts  on  the  lateral  surfaces  of 


'  Due  to  the  synunetry  and  uniformity  of  the  fields,  only  the  r  dependence  is  explicitly  indicated.  The 
superscript  c  denotes  the  confining  cylinder. 
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the  cylinder  regardless  of  whether  the  (a)  is  greater  than  zero  or  not  and  hence,  does 
not  prevent  the  barreling  instability  in  the  compression  specimen.  On  the  other  hand,  the 
confined  compression  configuration  that  is  suggested  here  ensures  that  the  cylindrical 
deformation  is  stable.  The  constraint  imposed  by  the  confining  cylinder  is  characterized 

by  the  quantities  cr^^(a)  and  u^(a),  and  is  fully  determined  by  the  geometry  and 

material  properties  of  the  cylinder.  If  the  hoop  strain,  SQQ^h)  =  at  the  outer  surface  of 

the  confining  steel  cylinder  is  measured  using  a  strain  gauge,  then,  the  radial  component 
of  displacement  and  stress  at  the  inner  surface  of  the  confining  cylinder  can  be  calculated 
in  terms  of  the  measured  hoop  strain  from  the  Lame  solution  (see  Timoshenko  and 
Goodier,  1934): 


{l-v,)a  +  {\  + 


(5) 

(6) 


and  are  the  modulus  of  elasticity  and  Poisson’s  ratio  of  the  confining  cylinder 

respectively.  Since  the  confining  cylinder  is  free  to  deform  in  the  axial  direction,  is 
taken  to  be  zero  and  the  Lame  solution  for  plane  stress  is  used  in  the  above  equations. 
Note  that  the  above  evaluation  of  the  constraint  assumes  complete  elastic  deformation  of 
the  confining  cylinder,  which  must  be  ensured  through  proper  choice  of  material  and 
geometry  for  the  confining  cylinder.  This  is  not  strictly  necessary;  in  fact,  allowing  the 
confining  cylinder  to  become  plastic  will  provide  a  way  of  altering  the  confining 
pressure,  but  will  also  necessitate  a  new  confining  cylinder  for  each  specimen.  As 
discussed  above,  continuity  conditions  at  the  interface  between  the  confining  cylinder 
and  the  specimen  require  that  the  radial  component  of  displacement  and  stress  be  equal. 
Thus,  equations  (5)  and  (6)  are  also  the  corresponding  displacement  and  stress 
components  in  the  specimen  at  r  =  a.  Note  that  the  contact  is  assumed  to  be  frictionless; 
this  can  be  accomplished  by  suitable  surface  finish  on  the  contacting  surfaces,  through 
proper  lubrication  and  by  minimizing  the  length  of  the  specimen.  For  the  axisymmetric. 
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homogeneous  deformation  the  strain  displacement  relations  imply  that  (a)  =  (a) ; 

furthermore,  equilibrium  equations  require  that  the  radial  and  hoop  stress  components  be 
equal:  cr^ia)  -  aggia) .  Thus,  the  complete  state  of  stress  and  strain  in  the  cylindrical 
specimen  can  be  determined  from  the  experimental  measurements  of  the  axial  strain,  , 
the  axial  stress,  ,  and  the  hoop  strain  in  the  confining  cylinder,  S|^ .  The  state  of  strain 
and  stress  in  the  specimen  is  homogeneous  and  is  given  below  in  terms  of  these 
measured  quantities: 


s 

^rr  ~  ^00  “  '~7r  )  Y 

a 

^zz  ~  *  *  -  iniposed  by  the  loading  frame  and  measured 

^rr  ~  ^60  ~  2 

^zz  -^a'" measured  using  a  load  cell 


Thus,  in  this  experimental  configuration  the  measurements  provide  a  complete 
determination  of  the  stress  components  <t^  =  agg  and  a.,  as  well  as  the  corresponding 

strain  components  =  Sgg  and  £..  without  invoking  any  description  of  the  constitutive 

behavior  of  the  specimen  material!  Since  these  are  the  principal  components  of  stress  and 
strain,  these  measurements  enable  a  simple  characterization  of  the  constitutive  response 
of  the  material.  This  technique  was  developed  in  an  earlier  work  and  has  been  quite 
useful  in  characterizing  the  mechanical  response  of  polymers  under  multiaxial  loading 
(Ravi-Chandar  and  Ma,  2000) 

This  configuration  provides  a  confining  compression  for  the  polyetherether 
ketoone  (PEEK)  specimens.  By  altering  the  geometry  and  material  of  the  confining 
cylinder,  the  confinement  can  be  varied,  our  recent  experiments  have  shown  that  under 
confined  compression,  PEEK  can  sustain  a  very  high  shear  stress  level.  The  plot  below 
indicates  the  axial  stress  vs.  axial  strain  variation;  while  the  unconfined  PEEK  exhibits 
shear  yielding  at  about  25  ksi,  the  confined  specimen  (in  two  repeat  trials)  is  shown  to  be 
capable  of  sustaining  a  much  higher  level  of  compression.  The  enhancement  in  the 
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capacity  of  PEEK  depends  on  the  level  of  the  confinement;  in  the  example  shown  the 
confinement  was  unreasonably  high  (about  70  ksi),  but  realistic  levels  of  about  30  ksi  can 
be  obtained  by  suitably  altering  the  microstructure  of  the  composite.  We  are  working  on 
translating  this  idea  to  the  design  and  fabrication  of  fiber  composites  by  generating  such 
compression  during  processing.  This  is  anticipated  to  lead  to  a  much  higher  tensile 
strength  in  the  composite. 


Axial  Strain 


Figure  6.  Variation  of  tlie  axial  stress  with  a.xial  strain  for  PEEK  un  uniaxial  and  confined  compression 

tests. 
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